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INTRODUCTION 

Let  X(t) ,  t  ^  0,  be  the  Ornstein-Uhlenbeck  process:  the 
stationary,  mean-zero  Gaussian  Markov  process  with  covariance 
Cov(X(t)  ,X(s) )  =  (|■)exp(-|t-s| )  .  Let 

Z(T)  =  si  x2(t)dt.  (1) 


In  this  paper  we  calculate  the  Laplace  transform  of  the 
probability  density  function  of  Z(T)  using  functional  integral 
methods,  invert  this  transform  by  expanding  it  in  an  infinite 
series,  and  study  the  asymptotic  behavior  of  the  density  and 
distribution  functions  as  T  -»  ».  The  functional  integral  approach 
makes  possible  a  parallel  analysis  of  the  density  and  distribution 
functions  of  Z(T)  conditioned  on  the  event  X(T)  =  0,  which  are  of 
interest  in  certain  applications  (Dankel,  to  appear) .  The  asymp¬ 
totic  limits  are  all  closely  related  to  the  inverse  Gaussian 
distribution  with  parameters  (T/2,T/2)  (Johnson  and  Kotz,  1970a) , 
whose  density  function  is,  to  within  a  constant  factor,  the 
asymptotic  form  of  the  conditional  density  of  Z (T) ,  but  not  of  its 
unconditional  density  (see  Theorem  3  below) . 

Z(T)  is  an  example  of  a  quadratic  form  in  noirmal  variables,  a 
subject  that  has  been  widely  studied;  see  (Johnson  and  Kotz, 
1970b)  for  a  comprehensive  survey  of  work  up  to  1970.  The  usual 
method  of  analysis  expresses  the  characteristic  function  of  the 
distribution  as  an  infinite  product  involving  the  eigenvalues  of 
the  covariance  operator.  Using  this  method  to  study  a  process  very 
similar  to  our  Z(T),  Stepian  (1958)  recognized  the  infinite 
product  as  an  entire  function  expressible  in  terms  of  elementary 
functions,  and  inverted  the  Fourier  transform  via  contour 
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integration  and  numerical  integration.  By  contrast,  our  functional 
integral  method  bypasses  calculation  of  eigenvalues  and  infinite 
products,  yielding  explicit  finite  formulas  for  transforms 
directly;  moreover,  the  terms  in  our  series  for  the  inverse 
transforms  (i.e.,  the  conditional  and  unconditional  probability 
distribution  functions  of  Z(T))  are  given  as  explicit  expressions 
involving  standard  higher  transcendental  functions.  The  explicit 
form  of  these  terms  makes  asymptotic  analysis  feasible. 

The  approximation  of  the  probability  density  function  of  Z(T) 
by  an  inverse  Gaussian  density  was  heuristically  introduced  by 
Grenander,  Poliak,  and  Slepian  (1959) ;  however,  they  gave  no 
definite  conditions  for  the  validity  of  this  approximation.  In 
Theorem  3  below,  we  show  that  the  approximation  is  asymptotically 
exact  to  within  a  constant  factor,  but  only  if  one  conditions  on 
the  event  X(T)  =  0.  As  far  as  we  are  aware,  the  study  of  such 
conditional  distributions  of  quadratic  forms  in  normal  variables 
has  been  a  field  of  research  largely  unexplored. 

We  shall  require  a  number  of  formulas  for  Laplace  transforms 
and  special  functions;  these  are  given,  along  with  references,  in 
the  appendices. 
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TRANSFORMS 

In  the  sequel  F^(t)  will  denote  the  probability  distribution 
function  of  Z (T) ,  and  f^(t)  the  corresponding  probability  density 
function.  We  shall  use  the  superscript  "c"  to  denote  quantities 
conditional  on  the  event  X(T)  =  0;  thus,  F^(t)  and  f^(t)  will 
denote  the  conditional  distribution  and  density  functions  of  Z (T) , 
respectively.  We  shall  denote  the  Laplace  transform  of  a  function 
f(t)  ,  t  €  [0,«)  ,  by  i£(f(t))  (s)  . 


Theorem  1; 

(A)  if(f,j,(t))  (s) 

=  e'^/^[t-[(l+2s)2+(i+2s)"2jsinh[^T(l+2s)2j  +  cosh  j^T(l+2s)  ^  j  J  "2 

(B)  if(f°(t))(s)  =  e*^/^  ^(l+2s)"2sinh[T(l+2s)2|  +  cosh[^T(l+2s)2j  j"2 , 
for  all  complex  s  such  that  Re(s)  I  (-!•)  . 


Proof :  Define  w(t)  by 

w(t)  =  y/2  [(t+l)2  x(t-  log(t+l))  -  X(0) 


,  t  ^  0. 

Then  (Doob,  1942)  w(t)  is  distributed  as  standard  Brownian  motion. 
Solving  this  equation  for  log(t+l)),  substituting  into  the 
integral  in  the  equation 


if(f^(t))(s)  =  E[exp(-s  jjx2(r)dr)] 

=  E[exp(-  I  log(t+l))-(||Yy 


]• 


2T 

where  S  =  e  -  1,  we  find 


=  E[exp[-  f  <>^])- 


(2) 


S  =  e^T-i 
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We  now  make  use  of  the  following  theorem  of  Cameron  and  Martin 
(1945):  let  y  be  a  real  number,  and  let  p(t)  be  a  positive, 

continuous  function  on  [0,1].  Let  be  the  least  characteristic 
value  of  the  Sturm-Lionville  problem 

f  •  (t)  +  Xp(t)f(t)  =  0  (3) 

f  (0)  =  f  (1)  =  0.  (4) 

Suppose  X  <  Xq,  and  is  any  solution  to  (3)  satisfying 

f;^(l)  =  0.  (5) 

Then 


:|exp[|  J'QP(t)(w(t)  +  V2  y)^dtjj  =  - 


■fx(0) 


(The  factor  of  2  arises  from  the  now  non-standard  normalization 
E(w^(t))  =  t/2  used  by  Cameron  and  Martin.  Also,  the  argument  of 
the  exponential  on  the  right-hand  side  of  (6)  is  given  by  them  as 

y2[xxjp(t)dt  +  J-^p(s)f^(s)ds]  dt]. 

However,  using  (3)  and  (5)  it  is  not  hard  to  show  that  this 

2 

expression  is  ^ust  y  fj^  (0)/fj^  (0)  ,  as  claimed  in  (6).) 

Since  we  desire  to  use  (6)  to  evaluate  the  integral  in  (2) , 
which  has  upper  limit  S,  we  also  need  the  following  scaling 
result:  let  X,y  be  real  numbers,  let  S,a  be  positive  real  numbers, 
and  let  p(t)  be  a  positive,  continuous  function  of  t  €  [O,®) .  Put 

I(a,X,S,y)  =  E  exp[|  XoP(at)  (w(t)  +  V2  y)  ^dt]  ]  .  (7) 


Then 


I(a,X,S,y)  =  I(aS,XS^,l,y/Vs) . 


(8) 
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To  see  this,  make  the  change  of  variables  t  =  rS  on  the  right-hand 
side  of  (7)  to  get 


I(a,X,S,y) 


exp 


J*Jp(aSr)  + 


V2  v'|2 
>/S  ■> 


Since  w(rS)/>/s  is  distributed  as  w(r)  ,  (8)  follows. 

Returning  now  to  (2),  we  see  that  putting  p(at)  =  (at+l) 
enables  us  to  write  it  as 


(s)  = 


-3-  e-y'/2dy 
V2? 


I(l,-s/2,S,y/v/2)  e"y  /^dy 

v'Jjr 


[s. 


SS‘ 


.  1  -iL.1 
'  'V2sJ 


1 

VStt 


the  last  equality  arising  from  (8) .  But,  according 
assuming  s  is  real  and  nonnegative. 


(9) 

to  (6), 


I 


exp 


fg' (0}  zi 

lg(0)  2S 


]- 


where  g(t)  is  any  solution  to 

g' • (t)  -  (|  S^) (St+l)'^g(t)  =0,  t  €  [0,1] 


satisfying 

g' (1)  =  0. 

The  general  solution  to  (11)  is 

g(t)  =  c^(tS+l)2+^  +  C2(tS+l)2-^, 

where 

u  =  f  . 


(10) 

(11) 

(12) 

(13) 

(14) 
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Since  no  function  of  the  form  (13)  satisfies  the  boundary 
conditions  (4)  (because  doing  so  would  imply  (S+1)^^  = 
(l-2u)/ (l+2u)  ,  for  \JL  1  f)  ,  the  hypothesis  X  <  of  the  Cameron- 
Martin  theorem  is  satisfied.  The  boundary  condition  (5)  requires 


Using  (14) ,  (17) ,  and  (18)  in  this  formula,  together  with  some 

straightforward  algebraic  manipulation,  yields  part  (A)  of  the 
theorem. 

For  the  proof  of  part  (B) ,  we  note  that  the  time-reversed 
process  X(T-t)  is  again  an  Ornstein-Uhlenbeck  process,  and  that 

Z(T)  =  jjx^(t)dt  =  jJx^(T-t)dt. 
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These  observations  imply  that 

E[exp(-s  J*Qx2(t)dt)  lx(T)=y]  »  E[exp(-s  Xjx^  (t) dt)  1  X(0)  =y 

so,  by  (2),  (8),  (10)  and  (16) 

‘£(f^(t|X(T)=y) )  (s)  =  exp(2(B+^)y^)  . 

Putting  y  =  0  yields 

and  (17) ,  (14) ,  and  simple  algebra  lead  to  part  (B)  of  the 
theorem . 

While  the  above  argument  shows  the  equalities  (A)  and  (B) 
only  for  s  0,  the  analyticity  of  both  sides  of  these  equations 
establishes  their  equality  for  Res  ^  (“^)  r  completing  the  proof. 

□ 


INVERSION 

We  invert  the  Laplace  transforms  of  theorem  l;  our  results 
stated  in  terms  of  distribution  functions,  rather  than 
densities. 


Theorem  2:  For  all  t,T  >  0, 


(A)  P„(t)  =  (2a3)=  j  a„J„(T,t), 


n=0 


^  (2n-l) ! ! (2n) ! 


n 


2%! 


(19) 


where 


(20) 
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Jn(T,t)  =  x"[v  +  T(2n+t)3  exp(-  +  T(2n+t)  )^)G^(v)dv, 

and 

Gn(v)  =  i  X^cosh(r)  [Df2n_i(i'/2Tv^)  -  °-2n-l 
^_2n-l  ®  parabolic  cylinder  function. 

(B)  F°(t)  =  Tr"V^/2gt(T-t)  2  b  jC(T,t), 

n=0 

where  =  l/2^n!  , 

n 

Jn(T,t)  =  Xq[v  +  T(2n+t)]  exp(-  +  T (2n+t) )  ^ ) G^(v)  dv, 

and 

G°(V)  = 

V  v-r  . 

/  cosh(r)  X  — 

r*0  u=0  >/Tru(v-r“U) 


{D2n(2Vv-r-u)  +  D2j^(-2v'v-r“U)| 


Proof:  We  begin  with  the  proof  of  (A).  Let  fj^(T,t) 
function  such  that 


i£(f^{T,t))  (s) 


sinh(Ts) 


+  cosh(Ts) 


and  let  Fj^(T,t)  =  XQf^(T,r)dr,  so  that 

i£(F^(T,t))  (s)  S£(f^(T,t))  . 

Then,  by  algebraic  manipulation, 


i£{Fj^(T,t))  (s) 


g-Ts/2 

v^(s+l) 


(21) 

(22) 

(23) 

(24) 

(25) 

(26) 

dudr. 

be  a 

(27) 
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^-Ts  ^  ^  s  >  0,  for  such  s  we  may  expand  the 
(s+i)  I 


Since 

quantity  in  curly  brackets  in  a  power  series  to  obtain 


00  3 

a:(F^(T,t))  (s)  =  -(2n) 


-Ts ( 2n+t)  2(s-l)^^  ] 

>/s(s+l)^^‘*’^-^  ' 


(28) 


with  the  given  by  (20)  (we  use  the  convention  (-1)!!  =  1)  . 
Taking  the  inverse  Laplace  transform  term-by-term  and  using  (Al) 
and  (A2) ,  we  have  formally  that 


F  (T,t)  =  i  e'^/2  2  a  x  (t) 

n=0  [T(2n+t)/“) 


{°-2n-l^^^^t-T(4n+l)) 

I  ^ 

“  D!2n-i(“^^2t-T(4n+l))}.  (29) 

To  justify  (29)  rigorously,  let  z  =  iy,  y  real,  in  (A3),  obtaining 

n  rivl  =  ^  ~  r“  e"^y^e"^^/^x2"dx 

-2n-l^^y^  (2n)!  0  ®  ®  x  ax, 

so  that 


or 


(30) 


Hence,  the  absolute  value  of  the  nth  term  of  the  series  in 
(29)  is  bounded  by 

2 

®xP(t-T(2n+t)), 

Which  is  easily  seen  to  be  e^  times  the  nth  term  of  a  convergent 
series  (e.g.,  by  the  ratio  test).  Thus,  by  dominated  convergence, 
for  s  >  1  we  may  take  the  Laplace  transform  of  the  series  in  (29) 
term-by-term,  arriving  at  the  right  side  of  (28) .  Since,  by  (28) , 
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F^(T,t)  and  the  series  in  (29)  have  the  same  Laplace  transform  for 
s  >  1,  they  are  the  same  function,  as  asserted  in  (29) . 

By  Theorem  1  and  (27) , 

S£(F  (t))(s)  =  I  ie(f  (t))  (s)  =  f[-|— la:(f  (T,t))(s)  ^ _ 

^  Si  J-  JsW^s+T 

=  \\-¥-]  ^(F,  (T,t))  (s)]  ^ _  , 

LLs^-l-J  ^  s-W'as+l 

where  the  notation  "s  -»  V2s+1''  indicates  that  s  is  to  be  replaced 
by  V2S+1.  Using  (A4)  and  (A5) ,  this  becomes 

if(F^(t))(s)  =  5£(2  J’^cosh(u)F3^(T,t-u)du)  (VJiTT)  . 

Applying  (A6)  and  (A7)  to  the  right-hand  side,  we  get 

S£(F^(t))  (s) 

=  Sf  ^(2/7r)2t“^/2e“'‘^/^  ve”’^  J*QCOsh(u)  F^  (T,  v-u)  dudv  (s)  . 


Using  (29) ,  and  integrating  term-by-term  (which  may  be  justified 
by  an  argument  similar  to  the  one  justifying  (29)),  we  see  that 
the  inner  integral  in  the  above  formula  is 


,J!/2  «  v-T(2n+i-)  r  2  , _ 

is -  2  a  V  (V)  r  cosh(u)<D^_^  .  (i>/2  (v-u) -T(4n+1) ) 

'  n=0  "  tT(2n+t),»)  J  I  -2n  1 


-  ofan-i (‘1^2 (v-u) -T(4n+1) ) }du. 


SO  that 
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‘£(F  (t))(s)  =  5£[(2A^)2t"^/2e2(T-t)  j  j- 

^  n=0  T(2n+t) 


-vV2t 


v-T(2n+t)  /-  _ _ 

X  i  J*  cosh(u)|D_2j^_2^(i''^2  (v-u) -T  (4n+l) ) 

-  Df2n_i("^'^2  (v-u)  -T(4n+1)  )|dudv  , 


the  termwise  integration  justified,  as  usual,  by  (30) .  Replacing  v 
by  v+T(2n+2‘)  in  the  integral,  and  invoking  the  uniqueness  theorem 
for  Laplace  transforms,  we  complete  the  proof  of  part  (A)  of  the 
theorem . 

The  proof  of  (B)  proceeds,  mutatis  mutandis,  as  above.  Let 
f^(T,t)  be  a  function  such  that 


5£(f®(T,t))  =  e*^/^ 


sinh(Ts)  +  cosh(Ts) 


2 

r 


and  let  F°(T,t)  =  J^f J(T, r) dr . 

Then,  arguing  as  before  (28) ,  we  find 


s£(Fj(T,t))  (s)  =  2^e^/2 


Then  we  have  formally 


;  viiL  r 

n=0  (2n)  I  '• 


^^-Ts(2n+t)  (s-1)^  1 ^ 

‘  Vs  (s+1) 


-  o2^T/2  *  Vlii 


F:r(T,t)  =  2^0"' ‘  2  - 

n=0  (2n) ! 

But  by  (A5)  and  (A8)-(A10), 


"  e 


-Ts(2n+t-)  (s- 


n!L  1 


=  "7=^  kn< 

'‘Vs(s+1)  0  Virr(t-r)  ^ 


(t-r) 


Vs(s+1) 


2Vt^)  +  D2j^(-2Vt^)  jdr 


X  (-l)"/((2n-l)  !  !  20, 
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SO,  using  (A2) ,  we  see  that  (32)  is 

F°(T,t)  (33) 

.4.  m/2  “  .t-T(2n+^)  E  (t-T(2n+^)-r) 

=  (2ir)  ^6^/2  2  by  (t)  {  S  -  - 

n=0  "  [T(2n+i-) ,«) '^  0  Virr(t-T(2n+i-)-r 

where,  for  brevity,  we  have  introduced  the  notation 

(34) 

To  justify  (33)  rigorously,  we  have  from  (All)  that,  for  z  real, 

|D2„(Z)  I  ^  -r-'  2=’'+"  x="a-2=‘^x 

=  22n+t  ^2^4  .  j-n-t  ySne-y^dy, 

which  is  simply 

|D2n(2)l  i  (2n-l)i!  e^  ,  z  real.  (35) 

Using  this  estimate,  and  the  elementary  result 

fn  "■  dr  =  IT,  L  >  0,  (36) 

Vr (L-r) 

we  see  that  the  absolute  value  of  the  nth  term  of  the  series  on 
the  right  side  of  (33)  is  bounded  above  by 

bn  *  >/ir  •  2  •  (2n-l)l!  exp(t-T(2n+2') )  , 

which  is  e^  times  the  general  term  of  a  positive  convergent 
series.  Hence,  for  s  >  1,  we  may  take  the  Laplace  transform  of  the 
right  side  of  (33)  termwise,  so  (33)  may  be  proved  as  was  (29) . 
Using  (35)  to  justify  further  termwise  integrations,  we  complete 
the  proof  of  part  (B)  of  the  theorem  as  in  the  proof  of  part  (A) . 
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ASYMPTOTICS 

In  this  section  we  study  the  asymptotic  behavior  as  T  "  of 
the  conditional  and  unconditional  distribution  and  density  func¬ 
tions  of  Z (T) .  Our  main  tool  is  the  following  special  case  of  the 
extended  Laplace  method  (briefly,  ELM  below)  for  the  asymptotic 
evaluation  of  Laplace  integrals  (Erdelyi,  1956,  p.  37) .  Let  g(v) 
be  a  function  defined  on  (0,")  such  that  g(v)  bv^~^  as  v  ^  0, 
X  >  0,  and  such  that,  for  k  >  0 

J*  g(v)e  ^^dv 
0 

exists  for  sufficiently  large  positive  T.  Then 
00 

S  g(v)e“^'^dv  br(X)(kT)"^,  T  -»  ». 

0 

Theorem  3;  Let  fQ(T,t)  be  the  probability  density  function  of 
an  inverse  Gaussian  distribution  with  parameters  T/2,  T/2: 

fQ(T,t)  =  (8Tr)‘^Tt“^/^  exp[-  |  +  |  “  |t]  '  ^  >  0,  (37) 

and  let  FQ(T,t)  be  the  corresponding  probability  distribution 
function: 

FQ(T,t)  (38) 

=  (t)  [Erfc^-(t/2)2  +  T/(2(2t)2)|  +  e'^Erfc[^(t/2)2  +  T/ ( 2  ( 2t)  ^  )  j  j  , 

2 

where  Erfc(x)  =  1  -  Erf(x)  =  2ir  e  ^  dt. 

X 

Then,  for  fixed  t  >  0,  as  T 

(A)  F^(t)  ^  (8t/T)2FQ(T,t)  ,  8(ir)"2tT"^/2  exp[-  |  |  r  (39) 
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(B)  f^(t)  ^  (8t/T)^fQ(T,t)  (40) 

(C)  F^(t)  ^  22FQ(T,t)  ^  4Tr“2t2T"^exp[-  |  +  |  "  ft]  (41) 

(D)  fj(t)  ^  2^f^(T,t).  (42) 

Note  that  the  final  asymptotic  equivalences  in  (A)  and  (C)  follow 

from 

Fo(T,t)  ^  (8A)2t2T"^exp[-  |  +  |  -  |ij  ,  (43) 


which  is  a  consequence  of  (A13) .  The  plan  of  the  proof  of  Theorem 
3  is  to  study  via  ELM  the  asymptotic  behavior  of  Jj^(T,t)/ 
d/dt(Jj^(T,t) )  /  Jn^T,t)/  and  d/dt (J^(T, t) )  as  T  -*  ®  (see  Lemma  2 
below),  which  is  turn  depends  on  the  asymptotic  behavior  of  G^(v) 
and  G®(v)  as  V  -»  0,  established  in  Lemma  1  below.  Then,  using  the 
uniformity  estimates  of  Lemma  3,  we  see  that  the  series  from 
Theorem  2  for  F^(t)  ,  f^(t)  ,  F^(t),  and  f!^(t)  are  all  asymptotic  to 
their  first  terms,  and  the  proof  is  completed  by  another 
application  of  the  n  =  0  case  of  Lemma  2. 


(A) 

(B) 


Lemma  1;  For  all  n  ^  0, 
Gn(v)  .  v^/^,  V  0, 


where  c^  =  (2n-l)  !  !nJ/3  ( (2n)  ! )  ^ . 

G^(v)  ^  d^v,  V  ^  0, 

where  d^^  =  2ir2  (-i) ” (2n-l)  ! !  . 


(44) 

(45) 

(46) 


(47) 
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proof ; 

(A)  Since  i  d/dy[Df2n-i(iy)  -  ^ °-2n-l 

as  y  -»  0,  it  follows  from  (22)  that 

Gn(v)  .  (-2'^^V3)D_2j^_^(0)D:2n_i(0)v^/2,  V  -»  0.  (48) 

But  it  is  an  easy  consequence  of  (A3)  that 

D_2n_i(0)  =  (ir/2)2(2n-l)  !!/(2n)  !, 
and  from  (A3)  and  (A12)  it  follows  that 
D:2n_i(0)  =  -2%!/(2n)l  , 

so  that  (44)  and  (45)  follow  from  (48) . 

(B)  By  (36)  ,  as  v-r  -»  0, 

f  - - -  Id _ (2Vv-r-u)  +  D^_  (-2>/v-r-u)  Idu 

u=0  VTU(V-r-u)  I-  > 

— •  ir®-2D,  (0)  =  (-l)"(2n-l)  !  !  , 

the  last  equality  resulting  from  (All).  Equations  (46)“(47) 

then  follow  from  (26) . 

□ 


Lemma  2;  For  fixed  t  >  0  and  for  all  n  ^  0,  as  T 


(A) 

Jn(T,t) 

^  c^r  (5/2)  ^2n+t)  V2t  (2n+i-) ) 

(49) 

(B) 

d/dt(J^ 

(T,t))  ^  2"^c^r(5/2)t2e"'^  (2n+t)  /2t(,p(2n+i-)  )2 

(50) 

(C) 

J^(T,t) 

dj^t^e"*^^  )  V2t^^^2n+t) )  "^ 

(51) 

(D) 

d/dt(J^ 

(T,t))  ^  2"^d^e"'^^^2”‘^^^^/^^(T(2n+t))  , 

(52) 

where  c_  and  d_  are  as  in  Lemma  1. 
n  n 


Proof !  From  (21)  we  have 


dl+Ij), 
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where 

and 

I2  =  T(2n+t)  J'Qe"’’^^”'^^^^/^e“'^^/^^G^(v)dv. 

By  ELM  and  Lemma  1(A)  ,  as  T  -♦  », 

I2  .  Cj^r{5/2)t^/2(T(2n+t))”^/^ 
and  =  0(t“^/^) ; 

thus  +  I2  ^  12#  and  part  (A)  follows. 

To  prove  (B) ,  differentiate  (21)  with  respect  to  t  to  obtain 
d/dt(Jn(T,t) )  =  (2t^)"VQ[v+T(2n+t)  ]^exp(-(v+T(2n+i-) )  V2t)G^(v)dv, 
the  differentiation  under  the  integral  sign  being  justified  by  the 
uniform  convergence  of  the  above  integral  for  t  in  any  compact 
subset  of  (0,«) .  Rewriting  the  above  expression,  we  have 


d/dt(J^(T,t)) 

Expanding  the  cube  under  the  integral  sign  and  investigating  the 
asymptotic  behavior  of  the  four  terms  as  T  -»  ®  by  ELM,  we  see,  as 
in  the  proof  of  (A)  above,  that  the  term  of  the  cube  with  the 
highest  power  of  T  dominates  the  others.  Hence, 


d/dt(J^(T,t)) 

(T(2n+t))^(2t^)”^e"'^  (2n+t)  /2tj.®g-T(2n+t) v/t^-v 
(T(2n+t))^(2t^)"^  e""^  (2n+t)  /2t  c^r  (5/2)  (T(2n+t) /t) 

=  2“^Cj^r(5/2)t2  (2n+t)V2t  (T(2n+t))^, 

which  is  (B) .  The  proofs  of  (C)  and  (D)  are  similar  to  those  of 
(A)  and  (B) ,  using  (25)  in  place  of  (21)  and  Lemma  1(B)  instead  of 


Lemma  1  (A)  . 
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Lemma  3 ;  There  exist  positive  sequences  i  =  1/2,  such 


that 


2  converges. 


n=0 

CO 


2  k^b^  converges. 


n=0 


(53) 


(54) 


and  such  that  for  fixed  t  >  0,  for  sufficiently  large  positive  T, 
we  have  for  all  n  ^  0: 


(A) 


(B) 


Jn(T/t) 


Jo(T,t) 


d/dt(J^(T,t)) I  ^  ^ 

d/dt(JQ(T,t)) 1  ^  ^n  ' 


(55) 


(56) 


(C) 


Jn(T/t) 


Jo(T/t) 


^  k2 


(57) 


(D) 


d/dt(j;;(T,t)) 


d/dt(JX(T,t)) 


^  k: 


n 


(58) 


Proof:  To  prove  (A),  we  may  begin  with  (21),  square  the 
binomial  in  the  argument  of  the  exponential,  and  drop  the  cross 
term  to  get 

|Jn(T,t)l  ^  /q[v  +  T(2n+t)]exp(-(v^  +  T(2n+i-)  ^)/2t)  |Gj^(v)  |dv. 

But  from  (22)  and  (30)  it  follows  that 

2 


so  we  have 


|Jn(T,t)|  i  +  T(2n+t)l2], 


(59) 


with 
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By  Lemma  2(A),  for  sufficiently  large  T,  we  have 

|jQ(T,t)  I  > 

where  here,  as  we  will  do  for  the  remainder  of  the  proof  of  Lemma 
3,  we  denote  by  K  any  expression  which  may  depend  on  t,  but  not  on 
T  or  n.  Hence,  for  such  T,  we  have,  for  all  n  ^  0, 


i  Ke^^T^/^e"’’  /2t(4n  +2n)  ^  T(2n+t)l2)- 

.  2 

Provided  T  /2t  >  l,  this  yields 

^  Ke^e“'*"  +  T(2n+i-)  I2) ]  • 

But,  if  T/2t  >  5/4,  T/2t  >  1  +  l/4n,  n  ^  1,  so  T^n/t  >  T(2n+i-)  , 
n  ^  1.  Thus 


i  Ke^e"'^”  (I^  +  T(2n+t)  I2)  ]  ,  n  ^  1. 

•  »  “X 

For  sufficiently  large  x,  the  function  xe  is  decreasing,  so  for 
sufficiently  large  T,  we  have  for  all  n  ^  1 


Jn(T,t) 

Jo(T,t) 


1  Kej^e"^^^[T^/^e"’^/2(I^  +  (T/2)l2)]. 

Since  the  last  factor  above  is  <  1  for  sufficiently  large  T,  part 
(A)  of  the  Lemma  is  proved:  putting 

2 

kj  =  1  and  =  Ke^^e"'^’’  ,  n  1, 

CO 

we  easily  see  that  2  kjljla  |  converges.  The  proof  of  part  (B)  is 

n=0  " 

very  similar.  To  prove  (C)  and  (D) ,  use  (26),  (35),  and  (36)  to 
derive  the  estimate  |g^(v)  |  i  27r*  (2n-l)  ! !  ve^  s  replace 
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(59)  and  argue  as  above,  putting 


, 2  -4n 

=  Ke„e 
n  n 


n  ^  1, 


and  checking  that  2  k"  b^l  converges.  The  proof  of  Leinina  3  is 

n=0  ” 

complete .  □ 


Turning  now  to  the  proof  of  Theorem  3 ,  part  (A) ,  we  have  from 
Theorem  2(A)  that 

lim  e“'^/^F  (t)/J  (T,t)  =  (2/Tr^)2  lim  2  a. 

T-^ 


=  (2/ir-')2 


3v2 


n=o  ^ 

J„(T,t) 


n 


n:o  ' 


the  interchange  of  the  order  of  limit  and  summation  being  allowed 
by  dominated  convergence  and  Lemma  3(A).  But,  by  Lemma  2(A), 


lim 

T-»® 


Jo(T,t) 


0,  n  ^  1, 


so,  since  a^  =  1,  we  have 

F^(t)  ^  (2/ir^)2  t"^/2^-t/2^T/2j^^^  ^ 

Applying  Lemma  2(A)  again  for  n  =  0  proves  part  (A)  of  Theorem  3. 
To  prove  part  (B) ,  we  see  from  Theorem  2(A)  that 


_d 

dt 


00 

2 

n=0 


^n  dt 


(Jn(T,t))  , 


the  termwise  differentiation  justified  by  Lemma  3(B).  Then 


lim 

T-»« 


_d 

dt 


(2A^)-^  t3/2e-V2(T-t)^^^^j 


1/— 

J'dt 


(Jo(T,t)) 


T-Ko  n=0 
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“  ni*"  ii:  (Jo(T.t))]  =^0  =  1' 

where  we  have  used  part  (B)  of  Lemma  2  and  Lemma  3.  Hence 


But,  by  part  (A) of  the  theorem  and  Lemma  2(B)  for  n  =  0, 
F^(t)  -  o(e’'/2  gf  (Jo(T,t))], 


SO  we  have 


(2/ir^)'^  ^  e'^/2  ^  (j^(T,t)), 

which,  by  another  application  of  Lemma  2(B)  for  n  =  0,  is  just 
part  (B)  of  the  theorem.  The  proof  of  part  (C)  is  exactly  parallel 
to  that  of  part  (A) ,  and  the  proof  of  part  (D)  is  exactly  parallel 
to  that  of  part  (B) ;  we  omit  the  details. 

□ 
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APPENDIX 

In  this  section  we  list  the  formulas  for  Laplace  transforms 
and  special  functions  employed  in  the  text  and  give  references  to 
standard  works  where  they  may  be  found,  using,  in  the  style  of 
Gradshteyn  and  Ryzhik  (1980) ,  the  following  abbreviations: 

GR  Gradshteyn  and  Ryzhik  (1980) 

OB  Oberhettinger  and  Badii  (1973) 

Thus,  "OB  234  (3.65)”  refers  to  formula  number  3.65  on  page  234  of 
Oberhettinger  and  Badii  (1973) .  The  formulas  are  listed  in  the 
order  in  which  they  are  needed  in  the  text.  Always  in  this 
appendix,  g(s)  denotes  the  Laplace  transform  of  f(t): 

g(s)  =  5£(f(t))(s)  =  /“  e"®4(t)dt. 

(AD  . 

OB  234  (3.65) 

(A2)  i^[xfb/a.«)  f(at-b)]  =  a'V^®/®g(s/a),  a,b  >  0, 

where  (t)  denotes  the  characteristic  function  of  the 

[c,®) ' 


interval  [c,~) :  x(t)  =  {j'  t  J  c}  ’ 
OB  3  (1.7) 


-zV4 


.00  ^-zx-x  /2  ^-p-1, 
GR  1064  (9.241(2)) 


(A3)  Dp(z)  =  -  J*o  e  X  ^  -dx,  p  <  0, 


(A4)  (cosh(t) ) (s)  =  s(s^-l)“^,  Re(s)  >  1 
OB  84  (9.2) 
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(A5)  93^(3)92(3)  =  ^[Sl  f^(u)f2(t-u)du]  (3) 

OB  7  (1.31) 

(A6)  9(a3+b)  =  (t/a)  j  (3)  ,  a  >  0 

OB  207  (1.3) 

(A7)  9 (Vi)  =  ie[i-  Tr“2t“V2  J.» 

OB  210  (1.27) 


(A8) 


js-JLl 


n 


=  r  (2-n)ir 


n+2 


(s+a) 

OB  182  (17.76) 


•Ij-n-l  ^j’ 


t  M(D2„(2(at)»)  + 


Re(s)  > 


(A9)  r(t-n)  =  (-l)'’2”TrV(2n-l)  !! 

GR  938  (8.399(3)) 

(AlO)  =  S£[(irt)“2j^  Re (3)  >  0 

OB  16  (2.27) 

(All)  Dp(z)  =  ir"2  gzV4  p 

GR  1064  (9.241(1)) 

(A12)  D^(Z)  =  pDp_3^(z)  -  t-2Dp(z) 

GR  1066  (9.247(2) ) 

_x2 

(A13)  Erfc(x)  ^  V  /x,  x  -»  +® 


-2(at)2)]] 
-Re (a) 


GR  931  (8.254) 
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